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Abstract 

We study a generalization of Swendsen-Wang algorithm suited for Potts mod- 
els with next-next-neighborhood interactions. Using the embedding technique 
proposed by Wolff we test it on the Symanzik improved bidimensional non- 
linear a model. For some long range observables we find a little slowing down 
exponent (z ~ 0.3) that we interpret as an effect of the partial frustration of 
the induced spin model. 
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I. INTRODUCTION 



The Swendsen-Wang algorithm [3] is known to be a very efficient way of generating 
configurations for a Monte Carlo simulation, owing to its little slowing down exponent. 
For the bidimensional Ising model the numerical data are consistent with r oc £-°- 35 ? to be 
compared with r oc £- 2J of the conventional, local algorithms (see for example [4-8]). 

However, the original formulation in the framework of the Potts spin model can't be 
easily generalized to other statistical systems, for example lattice gauge theories, in spite of 
many efforts to do it [9,10]. 

Wolff [11-13] has shown that it is possible to incorporate the Swendsen-Wang dynamics 
in a O(N) invariant, multicomponent statistical system, embedding Ising variables in the 
continuous degrees of freedom. This method proved to be extraordinarily efficient, with an 
almost complete absence of critical slowing-down r oc £« 0-1 . 

Motivated by a concrete application to the non linear, Symanzik improved 0(3) bidi- 
mensional sigma model, we study a simple generalization of the algorithm. 



II. THE ALGORITHM 



The basic idea of the Swendsen-Wang procedure is to introduce some auxiliary degrees 
of freedom in the model one wants to simulate. We apply the algorithm to the sigma model 
defined in our case by the Symanzik tree level improved action 



n,n 



(1) 



and by the constraint ip a (n)ip a (n) = 1. We fix an arbitrary unit vector r and parameterize 
the field as 



(p(n) = (p±(n) + \<p(n) ■ r\a n with a n 



ip(n) • r 



±1. 



(2) 



\<p(n) • r 

At fixed <p± and \ip ■ r\ the system is equivalent to a disomogeneous Ising model 

s = Y, J K n ']<W„/ ( 3 ) 



with a ferromagnetic nearest neighbor coupling and an antiferromagnetic third neighbor one 

/ 4 1 \ 

J[n,ri] = 2P\<p(n) ■ r\\<p(ri) ■ r\ ^ {--5 n >,n+n„ + y^„ >+27V J . (4) 

The partition function can be written, neglecting an irrelevant multiplicative constant, as a 
product of terms 



Z = J2U( 1 + k (n,n')S w ) with *( 

{a} n,n' 



n, n') 



y -J[n,n'] _ ^ 



(5) 



In each of them we introduce a new degree of freedom l n y, which can get values in {1,0}, 
so that in the general case we can write (S ai b = 1 — S a ^) 
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l + k{n,n')S Un ^ n , = Y, ( w hi( n ^ n ') S ^,a n ,Si nn ,,i + Wo fi {n,n')5 Un!Un ,S lnn , fi 

+ W lfi S an , an ,S lnyfi + W ,iS an , an ,S lnnl!l ) (6) 

with the conditions + W\fl = 1 + k and PVo,o + Wo,i = 1. As the Wjj constants are 
proportional to probabilities they must be non-negative. 

The variables we have added are in a one-to-one correspondence with the interactions 
of the model, and we sum over all configurations {/, a}. We obtain a new partition function 
which describes the joint dynamics of all the degrees of freedom. 

We start considering the evolution of the set {a} at fixed {/}. If we make the choice 

W lfi {n,ri) = W 0fi {n,ri) (7) 

the interaction between the sites n and n' becomes irrelevant if l n y = 0. This means that 
the spin system decomposes in a set of independent clusters Cj, each of them made of all 
the lattice sites which can be joined by a chain of l n y = 1 interactions. 
Inside each C\ the dynamics is described by an Ising-like effective action 

c - i ( Wi.ifoX) ^ c , Q X 

= ~ „Bc g Uu («,»')-*(«,»') J K '*< <8) 

which can be simplified by imposing another condition. If k > (i.e. if the interaction is 
ferromagnetic) we can choose W\ t \ = k, obtaining a model in which the two spins must be 
aligned in order not to pay an infinite action tribute. If the action is antiferromagnetic the 
analogous choice is W\ t \ = 0: in this case the two spins must be necessarily unaligned. 

For a fixed {a} configuration the probability distribution for l n y depends only on the 
two spins a n and o n < . For a ferromagnetic interaction the relevant term is (cfr. equation (6)) 

$<t,<t> [kS t ,i + Si fi ] + 6 a , a 'Si fi . (9) 

It follows that if the two spins are unaligned I must be set to zero. In the case of align- 
ment there is on the contrary an "activation" probability proportional to k/(l + k). If the 
interaction is antiferromagnetic we obtain 

5 a ,o> [(1 + k)6 lfi - kS lA ] + (1 + k)S a ^S lfi . (10) 

In this case if the spins are aligned it follows necessarily I = 0, in the other case we have 
/ = 1 with probability —k. 

From these considerations it follows that after the generation of the {/} set inside each 
cluster the spins automatically satisfy the constraint imposed by the equation (8), and that 
the only possible moves are the flippings of a cluster as a whole. 

In conclusion we can sum up the procedure as follows. After choosing a random direction 
r we set the / values with the appropriate probabilities. Next we construct the clusters, and 
flip each of them with some assigned probability. 

In absence of Symanzik improvement there are only ferromagnetic couplings, so each 
cluster is composed of aligned spins. In our case it is possible for two or more clusters of 
this type with opposite spin orientation to be joined by an antiferromagnetic active I. 
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This fact has two interesting consequences. First of all at f3 = oo our algorithm is no more 
ergodic, as one can easily construct field configurations that are left unchanged by it, apart 
for a trivial global flip. To see this consider three spins in the sites n, 11 + 71^ and n + 2n^. If 
o« 7^ < J n+2n^ the antiferromagnetic bound l n , n +2n^ is surely activated, and the same must be 
true for one of the two ferromagnetic ones l n ,n+n^ or ln+n^,n+2n^, so that all the spins belong 
to the same cluster. If a n = cr n+nfi = a n+2nii both l n ,n+n„ and ln+ nfl ,n+2 nfl are activated and 
the three spins are connected again. Only if a n ^ c n+n(i ^ cr n +2n IJ , there is a probability 
that the spins can be changed independently, but this cannot occur for a sufficiently smooth 
configuration. So at j3 = oo the algorithm can change only high wavelength modes, while if 
the field configuration is smooth all spins are connected in one unique cluster, and the only 
possible update is a global parity. If (5 is big but finite we expect the formation of a large 
cluster which connects nearly all the sites, and then a reduced decorrelation. We emphasize 
that this is not the case for the non improved model, where also at j3 = oo the only stable 
configuration is that in which all the spins are aligned. 

Another point is that the mean size of the cluster is no more connected with the sus- 
ceptibility, as is the case without improvement where a Fortuin-Kasteleyn representation 
exists [14]. 

III. PERFORMANCES 

In order to test the efficiency of the generalized algorithm we have measured the inte- 
grated autocorrelation time for the observables [15] 



n,m 

F = 



U 2 = Y J 4> a {n)4> a {m) (11) 

n,m 

- (j2 e 27ri{ni - mi)/L (f> a (n)(f> a (m) \ + - fe e 2 * i( " 2_m2)/ V W'm) (12) 

2 \n,m J ^ \ njm J 

Ei = ^ + ^)- (13) 

From the mean values of M 2 and F one can easily evaluate the two point function at the two 
smaller momenta available on a finite lattice 

X = V- 1 <M 2 >=G(p)\ ipi=Q (14) 
x ' = V~ l < F >= G(p) . (15) 

These are "long distance" dynamical quantities (in particular x is the susceptibility) from 
which it is possible to calculate 



(x/x') - 1 

4sin 2 (7r/L) 



(16) 



which is a possible definition of correlation length in a finite volume. On the other side the 
mean value of Ei is connected to the short distance dynamics 
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E = V- 1 <E 1 >=G(n)| |B|=1 . (17) 

We have chosen the single cluster update scheme proposed by Wolff [11-13], and we have 
measured also the size of the flipped cluster iV c . In our parameter range the ratio £/L is 
always less than 0.5, and the asymptotic scaling regime is not yet reached. For example we 
observe at best a 15% discrepancy between our measured correlation length and the exact 
value predicted by the Bethe ansatz [16]. 

We report in Table I the integrated autocorrelation time for x,x' £i an d N c , extracted 
from a series of 10 6 consecutive cluster updates. In Table II we list the analogous results 
obtained using the over heat bath algorithm [17]. 

We have calculated the integrated correlations applying the self-consistent method pro- 
posed by Madras et al. [18], and we have checked the stability of the result. 

In order to evaluate the critical slowing down exponent for a given observable O we try 
to fit our data using the standard finite size scaling ansatz 

(18) 

Here £(/3, L) is the measured correlation length defined by (16) and <j> is a unknown universal 
function. As an example we report in Fig. 1 £~ z x versus £/L for all measures we have taken, 
using the value z = 0.3 which gives a reasonable result. 

Our best estimate for the critical slowing down exponents are reported in Table III 
and IV. As one can see the cluster algorithm performs certainly better in respect of the 
local one. For the long range quantities M 2 and F we argue that 0.2 < z < 0.4. It is 
interesting to note that for the local quantity E\ the results are consistent with a total 
elimination of slowing down. This is in some sense an intermediate situation between a local 
algorithm, which decorrelates short scales much better than long ones (see table IV for the 
over heat bath case), and the usual Swendsen Wang which reduces slowing down with the 
same efficiency at all scales. 

In Fig. 2 we plot the ratio between the measured susceptibility and the mean size of 
flipped cluster versus the correlation length. As we have anticipated with our generalized 
algorithm x is n °t proportional to N c , as one can easily see. We try to interpret the plot 
in the following way: for £ < 0.2 L the finite size effects are small, and we can see that 
the cluster size grows more rapidly than the "physical" size connected to the susceptibility. 
This is consistent with the discussion of the previous section about the expected behavior 
at large j3 values. For £ > 0.2 L volume effects prevent more effectively the cluster size than 
susceptibility from growing, hence N c /x decreases. 

IV. CONCLUSIONS 

Our results show that the proposed algorithm is effective in reducing the slowing down 
at short and long scales. In the last case the slowing down is not completely eliminated and 
we can interpret this fact in two equivalent ways. 

As there is not a proportionality between the cluster size and the physical scale of the 
model the algorithm is not forced to operate on the modes physically relevant. 



r = £(/?, L) z °c 
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From another point of view we have seen that the non optimal behavior at large f3 is 
connected to the simultaneous presence of ferromagnetic and antiferromagnetic interactions 
in the effective spin model, that becomes frustrated. It is well known that in presence of 
frustration the reduction of slowing down is an extremely difficult task. 

In our case the frustration is small, and the algorithm is in any case more efficient than a 
local one. We have worked out a more elaborate generalization of Swendsen Wang algorithm 
that could be effective in reducing the excessive growth of cluster size, and we are testing it 
to see if it is possible to further reduce slowing down in this model [19]. 

We are also extending our study to larger correlation length, in order to be sure that the 
dynamical exponents we have extimated are really the asymptotic ones. 
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FIGURES 



FIG. 1. Finite size scaling of r^ n * 
FIG. 2. The ratio N c /x versus the measured correlation length. 
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TABLES 



TABLE I. Results of numerical simulations for the autocorrelation times with Swendsen Wang 
algorithm. 
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P 


T int 

V 
A 




-int 
x' 




-int 
T Ri 


32 


1.250 


11.2 


(5) 


22.6 


(14) 


0.55 (6) 


32 


1.300 


12.1 


(6) 


23.9 


(16) 


0.51 (6) 


32 


1.350 


13.8 


(7) 


24.0 


(16) 


0.79 (1) 


32 


1.400 


14.6 


(7) 


24.6 


(16) 


0.50 (6) 


32 


1.450 


14.6 


(7) 


23.2 


(15) 


0.50 (6) 


32 


1.500 


15.9 


(8) 


21.1 


(13) 


0.50 (6) 


64 


1.250 


21.4 


(13) 


17.8 


(10) 


0.520 (6) 


64 


1.300 


14.3 


(7) 


17.3 


(10) 


0.500 (6) 


64 


1.350 


12.8 


(3) 


20.7 


(13) 


0.510 (6) 


64 


1.400 


13.4 


(7) 


26.2 


(18) 


0.500 (6) 


64 


1.450 


15.8 


(8) 


27.8 


(20) 


0.500 (6) 


64 


1.500 


16.7 


(9) 


25.5 


(17) 


0.500 (6) 


64 


1.525 


16.3 


(9) 


27.7 


(20) 


0.500 (6) 


64 
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17.4 


(10) 


26.1 


(18) 


0.500 (6) 


64 
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19.6 


(12) 


26.9 


(19) 


0.500 (6) 


64 


1.600 


17.6 
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25.3 


(17) 


0.500 (6) 


128 


1.300 


26.4 


(18) 


27.5 


(19) 


0.500 (6) 


128 


1.400 


17.5 


(10) 


19.9 


(12) 


0.500 (6) 


128 


1.500 


15.1 


(8) 


29.5 


(22) 


0.500 (6) 


128 


1.525 


15.6 
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30.1 


(22) 


0.500 (6) 


128 


1.550 


16.9 


(9) 


31.4 


(24) 


0.500 (6) 


128 


1.575 


17.0 


(9) 


31.8 


(24) 


0.500 (5) 


128 


1.600 


20.9 


(13) 


34.6 


(28) 


0.500 (5) 


128 
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(14) 
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(28) 


0.500 (5) 


128 
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22.0 


(14) 
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0.500 (6) 


128 
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(13) 


32.4 
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0.500 (6) 


128 
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(21) 


0.500 (5) 


128 


1.725 
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0.500 (5) 
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1.300 
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(31) 


0.500 (6) 


256 
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37.4 
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0.500 (6) 


256 


1.675 


22.2 


(14) 
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(35) 


0.500 (6) 


256 


1.700 


24.1 


(16) 


38.0 


(32) 


0.500 (6) 


256 


1.725 


25.3 


(17) 


44.5 


(40) 


0.500 (6) 


256 


1.750 


30.1 


(22) 


49.2 


(47) 


0.500 (6) 


256 


1.775 


27.4 


(19) 


39.6 


(34) 


0.500 (5) 


256 


1.800 


27.4 


(19) 


43.1 


(38) 


0.500 (6) 
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TABLE II. Results of numerical simulations for the autocorrelation times with over heat bath 
algorithm. 
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L 


P 


—int 

T x 


int 
V 


int 


32 


1.250 


8.8 (4) 


5.8 (2) 


2.37 (6 


32 


1.300 


10.9 (6) 


6.3 (2) 


2.70 (7 


32 


1.350 


1 O A ( T\ 

12.4 (7) 


T /I f 0\ 

7.4 (6) 


O TC / T 

2.75 (7 


on 

32 


1 /( a a 

1.400 


11.2 (6) 


TO / O \ 

7.3 (3) 


O AO ( f* 

2.48 (6 


32 


1 /irn 

1.450 
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9.8 (5) 


T O ( O \ 

7.6 (3) 


O 1 1 / f 

2.11 (5 


32 
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1.500 


O O / A \ 

8.8 (4) 


6.9 (3) 


1 AO ( A 

1.93 (4 


64 
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7.8 (3) 


5.9 (2) 


2.16 (5 
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1 OTA 

1.250 


11.1 (6) 
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7.7 (3) 


o o o ( o 

2.33 (6 


64 


1.300 


17.7 (12) 


10.1 (5) 


2.38 (6 


64 


1.350 


O A O / 1 A \ 

24.8 (19) 


1 A T /a\ 

14.7 (9) 


2.62 (7 


f* A 

64 


1.400 


O O A / O A \ 

33.4 (30) 


18.7 (13) 


2.60 (7 


f 1 A 

64 


1 1 7 A 

1.450 


07 o / n f\ 

37.8 (36) 
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23.6 (18) 
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2.58 (7 


O A 

64 


1 r" a A 

1.500 
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40.4 (40) 
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34.9 (32) 


OO c / 1 o\ 
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64 
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1.550 


OA 1 ( C\£} \ 

30.1 (26) 


OO 1 / 1 f* \ 

22.1 (16) 


o i A 1 r 

2.14 (5 


O A 

64 


1.575 


26.0 (21) 


OA A / 1 r \ 

20.9 (15) 


2.04 (5 


f 1 A 

64 


1 f* AA 

1.600 


OO O /O O ^ 

28.2 (23) 


O O 1 / 1 T\ 

23.1 (17) 


1 A f / /I 

1.95 (4 


128 


1.200 


7.0 (3) 


6.8 (3) 


2.38 (6 


1 O O 

128 


1 O A A 

1.300 


11.3 (6) 


1 A O ( O \ 

10.8 (6) 


O A f ( O 

2A& (6 


1 O O 

128 


1 A A A 

1.400 


1 f* O / 1 A\ 

16.8 (10) 


1 A f /A\ 

14.5 (9) 


O O O / T 

2.63 (7 


1 OO 

128 


1 r aa 

1.500 


AO T / r o \ 

48.7 (53) 


O T f* / O/^ \ 

37.6 (36) 


2.66 (7 


1 OO 
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1.525 


AO A /l O / \ \ 

93.0 (139) 
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121.3 (207) 


1 AO A / 1 T /I \ 
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1 AO / A 
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19.1 (13) 


19.1 (13) 


O A f ( O 

2A& (6 


o r o 

256 


1 A A A 

1.400 


£?0 O /T£?\ 

62.2 (76) 


/(TO /ri \ 

47.8 (51) 


n O A / T 

2.60 (7 


o r o 

256 


1 r a a 

1.500 


1 O /I t /o^r^ 

184.7 (275) 


AO O ( C\'l\ 

92.3 (97) 


o r o / r 

2.58 (5 


256 


1.600 


255.5 (447) 


183.6 (272) 


2.33 (4 


256 


1.625 


/(OO A / *l 1 O /I \ 

488.9 (1184) 


1 A P" A ( O A A \ 

195.4 (299) 


2.27 (4 


256 


1 f A 

1.650 


455.7 (1066) 


319.2 (625) 


O AO / O 

2.08 (3 


o r 

256 


1.675 


582.6 (1541) 


/I 1 O O /A1 0\ 

412.6 (918) 


o c\o ( o 

2.06 (3 


o r o 

256 


1 T A A 

1.700 


r T£? A /01 /I T\ 

576.9 (2147) 


oio 1 /or a^ 

313.1 (859) 


1 A£? / O 

1.96 (3 


256 


1.725 


373.6 (1119) 


247.2 (602) 


1.96 (3 


256 


1.775 


377.4 (1136) 


378.5 (1141) 


1.88 (4 


256 


1.800 


524.2 (1860) 


246.1 (593) 


1.76 (4 



11 



TABLE III. Critical slowing down exponents for the generalized algorithm. 



Li 


L 2 


Z M 2 


z F 


ZGi 


32 


64 


0.19 


0.11 


- 0.0 


32 


128 


0.24 


0.23 


- 0.0 


32 


256 


0.26 


0.25 


- 0.0 


32 


512 


0.27 


0.23 


- 0.0 


32 


1024 


0.32 


0.29 


- 0.0 


64 


128 


0.23 


0.30 


- 0.0 


64 


256 


0.24 


0.22 


- 0.0 


64 


512 


0.23 


0.18 


- 0.0 


64 


1024 


0.30 


0.26 


- 0.0 


128 


256 


0.25 


0.16 


- 0.0 


128 


512 


0.20 


0.10 


- 0.0 


128 


1024 


0.31 


0.21 


- 0.0 


256 


512 


0.24 


0.16 


— 0.0 


256 


1024 


0.35 


0.20 


- 0.0 


512 


1024 


0.47 


0.38 


- 0.0 




lABLh/ IV. Critical slowing 


; down exponents tor the over 


neat bath algorithm. 




Li 


L 2 




z F 




32 


64 


1.65 


1.63 


- 0.0 


32 


128 


1.65 


1.73 


- 0.0 


32 


256 


1.69 


1.76 


- 0.0 


04 


Izo 


1 Q 1 

l.oi 


l.yl 


^ u.u 


64 


256 


1.43 


1.83 


- 0.1 


128 


256 


1.66 


1.67 


- 0.2 



12 



10 

9 
8 
7 
6 

>< 
a 

^ 5 

d 



41 







0> 



n 



ii 



H 



H 



3 

2 

1 


0.15 




0.2 



0.25 



0.3 



0.35 



0.4 



0.45 



0.5 



